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SUMMARY 


The  specific  objective  of  this  program  was  to  establish  the  feasibility 
of  developing  an  analytical  model  for  solving  the  Reynolds  Equations  of 
Motion  for  fully  developed  two-dimensional  channel  turbulent  flow.  During 
the  course  of  these  studies,  feasibility  of  using  the  Turbulent  Flow  Models 
to  provide  the  closure  necessary  to  solve  the  Reynolds  Equation  has  been 
demonstrated. 

Closed  form  solutions  have  been  obtained,  by  use  of  the  Turbulent  Flow 
Model,  for  tJhe  mean  velocity  profile,  the  cross  correlation  (Reynolds)  stress 
term  and  for  the  root-mean-square  levels  of  the  axial  and  lateral  velocity 
perturbations.  These  solutions  compare  very  well  with  measurements  recorded 
by  John  Laufer  in  1951.  The  Turbulent  Flow  Model  was  also  used  to  obtain 
solutions  of  the  three  terms  of  the  Turbulent  Kinetic  Energy  Equation  - 
Production,  diffusion  and  dissipation  terms.  These  solutions  also  agree  well 
with  the  measured  data. 

The  studies  indicate  that  a  definite  variation  in  turbulence  eddy  properties 
exists  across  the  channel.  These  variations  can  be  determined  from  measurements 
of  the  flow  field  unsteady  pressure  fluctuations. 

A  procedure  for  extending  the  current  procedures  to  the  case  of  axisymmetric 
flow  is  included. 
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1.0  INTRODUCTION 


The  large  mass  flows  and  high  degree  of  integration  required  of  V/STOL 
type  propulsion  systems  affect  the  local  flow  fields  around  the  aircraft. 

These  propulsion  induced  effects  can  substantially  alter  the  performance  of 
the  aircraft  in  V/STOL  operations.  Forces  resulting  from  suckdown  and  jet 
fountain  can  modify  the  lift  and  pitch  and  roll  moments.  Such  forces  can 
be  appreciable  as  measured  in  scale  model  test  programs.  Similarly,  exhaust 
gas  reingestion  degrades  thrust  available.  Consequently,  it  is  important  to 
properly  model  the  propulsive  jet  when  testing  with  scale  models.  Different 
results  have  been  found  when  using  different  jet  sources  such  as  model  fan, 
compressed  air  jet,  or  ejector  type  propulsion  simulators.  Apparently, 
different  jets  which  are  similar  in  velocity  and  diameter  create  different 
downwash  and  ground  jet  flow  fields  with  differing  amounts  of  entrainment 
and  different  airframe  induced  forces  and  moments.  From  such  results,  past 
V/STOL  aircraft  performance  has  not  been  predicted  with  acceptable  accuracy 
even  though  extensive  controlled  model  data  was  used.  Improved  jet  entrain¬ 
ment  simulation  and  understanding  of  model  jet  scaling  is  required  to  over¬ 
come  this  type  of  V/STOL  program  hazard. 

Jet  entrainment  is  dependent  on  the  jet  velocity  profile  and  turbulence 
level  which  at  present  are  not  accounted  for  in  model  testing  or  in  mathematical 
models  of  jet  entrainment  and  mixing.  The  ultimate  goal  of  this  program  is  to 
develop  a  model  that  will  account  for  all  the  major  parameters  in  the  mixing 
process,  such  that  jet  entrainment  can  be  modeled. 

The  objective  of  this  program  was  to  establish  the  feasibility  of 
solving  the  Reynolds  Equations  analytically  through  application  of  the 
statistical  fluid  dynamic  model  of  turbulence  to  channel  flow.  This  fluid 
dynamic  model  has  been  successfully  applied  to  the  turbulent  flow  in  inlet 
systems  and  will  provide  the  link  between  the  statistical  properties  of 
turbulence  and  the  steady  state  velocity  gradients.  Results  will  enable 
prediction  of  the  steady  state  velocity  profile  and  Reynolds  stresses  through¬ 
out  the  channel  flow  field. 

The  technique  will  enable  an  analytical  solution  of  the  Reynolds  Equations 
of  Motion  without  recourse  to  empirical  constants  such  as  mixing  length  or  eddy 
viscosity  coefficient.  The  results  are  applicable  to  turbulent  flow  boundary 
layers  and  flow  mixing  problems  in  general,  and  ultimately  to  aircraft  propulsion 
induced  effects. 
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2.0  DISCUSSION 


rile  ultimate  goal  of  this  program  is  the  development  of  an  analytical  model 
for  the  simulation  of  jet  entrainment  and  mixing.  Jet  entrainment  can  be 
character ized  as  a  highly  complex  case  of  turbulent  free  shear  flow.  The  case 
of  fully  developed  turbulent  flow  in  a  two-dimensional  channel  is  the  simplest 


flow  problem,  and  therefore  was  chosen  as  the  model  for  the  present 
study.  A  previously  developed  turbulence  model'  was  used  to  obtain 
the  tuibulent  flow  equations  without  recourse  to  strictly  empirical 


feas  ib i 1  i  ty 
closure  of 
constants . 


The  model  has  physical  basis  in  the  coherent  vortex  approach. 


A  discussion  of  the  elements  investigated  during  the  feasibility  study  art- 
presented  in  this  section.  A  review  of  the  turbulent  flow  model  is  first  pre¬ 
sented.  This  model  Is  developed  using  fundamental  laws  of  fluid  mechanics  and 
siatistical  mathematics  and  is  based  on  a  series  of  eddies  of  random  size  and 
strength  convected  by  the  local  flow.  The  elementary  results  of  this  turbulence 
model  are  the  unsteady  velocity  perturbations.  These  perturbations  are  then 
combined  to  obtain  the  mean  and  root -mean-square  parameters  that  define  the  turui 
lent  flow  field.  The  variations  of  these  parameters  depend  on  the  statistical 
behavior  of  t he  turbulence  model  random  variables. 

These  random  variables  (eddy  size,  strength,  location,  spin  direction,  ei 
are  in  turn  dependent  on  the  actual  turbulent  flow  field.  The  second  part  of  Ur 
study  presented  in  this  section  shows  how  turbulence  flow  measurements  in  teiw> 
the  fluctuating  total  pressure  root-mean-square  values  define  the  distribution 
vortex  properties.  Since  these  measurements  were  not  available,  representat ' v 
distributions  were  used  and  refined  until  good  correlation  with  test  data  wa- 
oh ia i ned . 


Derivation  of  the  unsteady  velocity  correlation  terms  is  presented .  Ihas- 
terms  include  the  axial  and  lateral  velocity  components  (u'2  and  v'2)  and  ih 
cross  correlation  terms  (u 1 v 1 ) .  Also,  terms  in  the  Turbulent  Kinetic  Energy 
Euuation  are  presented. 

The  form  of  the  Reynolds  Equation  of  Motion  applicable  to  the  two-d  irnerrs  ion. 
channel  flow  is  presented.  A  discussion  of  the  boundary  conditions  and  (heir 
effects  is  also  included. 


Comparisons  are  made  between  the  various  flow  parameters  as  measurer*  by 
Lanier  and  reported  in  Reference  2  and  the  corresponding  predictions  made  by 
lb-'  current  analytical  work.  These  predictions  include  all  the  elements  of 
the  study. 

An  outline  of  the  steps  for  extending  the  analytical  method  to  the  case  c>! 
axi symmetric  flow  Is  presented.  This  i$  the  second  step  in  developing  the 
analytical  model  for  the  simulation  of  jet  entrainment  and  mixing. 


2.1  TURBULENT  FLOW  MODEL 


An  analytical  model  of  turbulent  inlet  flow  was  developed  and  reported 
in  Reference  1.  This  model  accounts  for  the  fundamental  fact  that  turbulence 
is  random  in  character,  and  consequently  that  a  statistical  analysis  of  the 
flow  is  an  essential  requirement.  To  model  this  random  inlet  turbulence,  it 
was  hypothesized  that  measured  pressure  fluctuations  result  from  a  random 
distribution  of  discrete  vortices  being  convected  downstream  by  the  mean  flow. 
Tangential  velocity  and  local  pressure  gradients  associated  with  each  vortex, 
when  superimposed  on  the  mean  flow,  produce  these  measured  pressure  fluctuations. 
Fundamental  fluid  mechanics  is  used  to  describe  the  characteristics  of  discrete 
vortices  and  their  superposition  on  the  mean  flow,  while  standard  statistical 
methods  are  used  to  generate  the  random  nature  of  the  flow  field.  By  use  of 
this  model,  the  size  and  strength  of  the  momentary  low  pressure  regions  are 
def i ned. 

A  physical  model  of  turbulent  inlet  flow  is  developed  by  consideration 
of  the  nature  of  turbulence.  This  model  is  then  formulated  mathematically  by 
use  of  fundamental  laws  governing  fluid  flow  and  by  application  of  statistical 
techniques  to  the  random  character  of  the  flow. 

Physical  Description  of  Turbulent  Flow  Model 

Fully  developed  turbulent  flow  in  a  channel  will  have  a  symmetrical 
velocity  profile  as  Indicated  in  Figure  1(a).  By  using  appropriate  flow 
visualization  techniques,  a  time  exposure  taken  of  this  flow  will  show  stream¬ 
lines  as  sketched  in  Figure  1(b).  Photographs  taken  by  Prandtl,^  with  the 
camera  moving  at  speeds  approaching  that  of  the  fluid,  show  vortices  being 
convected  along  with  the  local  mean  flow  as  in  Figure  1(c)  and  (d) .  There¬ 
fore,  the  concept  of  a  straight  streamline  is  only  true  on  a  time  averaged 
basis.  The  actual  fluid  particle  experiences  a  non-uniform  motion,  which  is 
dependent  on  the  turbulence  level. 

Roshko  discusses  the  case  of  the  turbulent  mixing  layer  in  Reference  A. 
Figure  2,  which  has  been  taken  from  Reference  4,  is  a  schematic  of  this  type 
of  flow  in  which  the  uniform  streams  on  either  side  of  the  mixing  layer  are 
characterized  by  their  velocities,  U]  and  U2,  etc.  Shadow  pictures  of  the 
flow  revealed  the  presence  of  well  defined  large  structures  (Figures  3  and  4) 
which  have  the  appearance  of  breaking  waves  or  rollers  or  vortices,  super¬ 
imposed  on  a  background  of  finer  scaled  turbulence. 

It  is  hypothesized  that  turbulent  channel  flow  is  composed  of  vortices 
of  various  size  and  strength  that  are  convected  downstream  at  the  local  velocity. 
This  hypothesis  is  fundamental  to  the  proposed  analysis  and  is  shown  schemat ica 1 1 y 
in  Figure  5.  Such  vortices  created  by  the  steady  state  velocity  profile  (U,  V) 
produce  unsteady  velocity  components  (u,  v)  as  measured  in  the  coordinate  system 
fixed  to  the  channel  walls.  The  flowfield  about  each  vortex,  however,  is  steady 
in  a  coordinate  system  fixed  to  the  vortex,  i.e.,  moving  at  the  local  velocity. 
Consequently,  the  vortex  flowfield  can  be  analyzed  using  steady  state  flow 
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c:  Streamlines  Vi  sable 
to  Observer  Moving 
About  1/3  U 


d:  Streamlines  Visable 
to  Observer  Moving 
at  U„ 


Figure  1.  Turbulent  Flow  In  a  Two-Dimensional  Open 
Channel  after  Photographs  by  Prandtl 


Figure  2.  Plane  Mixing  between  Two  Streams  with 
Velocities  and  Densities  and 
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Figure  3.  Mixing  Layer  Between  Helium 
(upper)  and  Nitrogen  Streams,  LL/U.  - 
O.38,  After  Roshko 


Figure  k.  Mixing  Layer  Between  Nitrogen 
(upper)  and  a  Helium-Argon  Mixture  of 
Same  Density,  After  Roshko 
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VELOCITY  PROFILE 


VORTEX  DISTRIBUTION 


Each  vortex  is  described  by  its  size,  a 

strength,  V0 

m 

location,  Y 
rotation,  +1 

Figure  5.  Vortex  Model  Flow  Field 


equations.  The  velocity  fluctuations  are  generated  from  superposition  of  the 
vortex  flowfieid  on  the  local  flow.  The  turbulent  nature  of  the  flow  results 
from  a  distribution  of  these  vortices  having  random  size,  strength,  location, 
and  orientation.  The  turbulence  terms  u'2,  v  ,  u'v'  can  be  determined  by 
application  of  statistical  methods  to  the  combined  flowfieid  definition. 

Isolated  Vortex  F 1 owf i e 1 d 


As  the  first  step  in  the  mathematical  development  of  the  flow  model,  the 
two-dimensional  incompressible  Navier-Stokes  Equations  of  Motion  were  solved 
to  derive  the  flow  properties  of  an  isolated  vortex  with  a  viscous  core  and 
inviscid,  potential  flow  at  large  radii.  Definition  of  the  variables  that 
comprise  the  flow  model  are  illustrated  in  Figure  6.  The  normalized  velocity 
field  associated  with  this  vortex  for  a  coordinate  system  located  at  the  center 
of  the  vortex  is  given  by  Equation  (1). 


/v 


On 


n (r/a) 


„-l/2[ (r/a;  -  1] 


(1) 


where  v  is  the  tangential  velocity  and  Ma"  the  radius  at  which  the  tangential 
velocity  is  a  maximum.  The  parameters  "a"  and  V-^  ,  are  used  to  characterize 
the  size  and  strength  of  the  vortex,  respectively,  with  "a"  being  considered 
the  radius  of  the  vortex  core.  The  associated  static  pressure  profile  through 
the  vortex  is  given  in  Equation  (2). 


P  -  P  =  -  £.  v  e  “t  (r/a)  -  1] 

o  2  0  me 

where  Pq  is  the  pressure  at  r/a  >>1 ,  and  p  is  the  density. 


(2) 


The  vortex  tangential  velocity,  Equation  (1),  can  be  expressed  in  terms 
of  the  horizontal  (axial)  and  vertical  (lateral)  velocity  components.  Since 
the  radius  to  a  point  (x,y)  in  the  flow  field  is  r^  =  x^  +  y2,  the  velocity 
components  are: 


u  =  -nv^  (y/a)e  'it(x/a)2  +  (y/a)2  -  1] 
v  =  nv(^  (x/a)e  "it(x/a)2  +  (y/a)2  “  1] 


(3a) 

(3b) 


Flowfieid  Due  to  Single  Vortex 

To  determine  the  effects  of  this  isolated  vortex  on  the  flowfieid  "fixed" 
to  the  channel,  a  coordinate  transformation  is  required.  Consider  the  sketch 
in  Table  i.  The  point  of  interest  in  the  flowfieid  Is  the  "total  pressure 
probe"  face.  It  is  desired  to  define  the  flowfieid  characteristics  at  the  probe 
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TABLE  I 


VELOCITY  AND  PRESSURE  PERTURBATIONS  DUE  TO  VORTEX  IN  2-D  DUCT 


II 


>i 

a 


-  Duct  Height 

-  Probe  location 

-  x  location  of  vortex  "i" 

-  y  location  of  vortex  "i" 

-  Vortex  size  (radius) 

-  Vortex  tangential  velocity 
at  radius  a 


v  -  Perturbation  velocity  in  y  direction 

u  -  Perturbation  velocity  in  x  direction 

U  -  Velocity  in  x  direction 

n  -  Index  of  rotation,  counterclockwi se=n  i 
clockwise  =  -1 


The  following  velocities  are  due  to  a  vortex  as  described  in  above 
sketch  and  were  developed  in  Reference  1. 


(1-1  > 


(1-2) 


(1-1) 


w 


due  to  the  presence  of  a  single  vortex.  At  any  instance  of  time  the  single 
vortex  is  located  at  some  distance  x. ,  y,  from  the  probe.  Since  the  vortex 
is  convected  by  the  local  flow,  its  location  in  the  x-direction  can  be  expressed 
as 


x.  -  Ut 

where  the  time  t  ■  0  is  taken  as  the  instance  when  the  vortex  passes  the 
probe  station  (x  ■=  0).  The  resultant  pressure  and  velocity  perturbation 
equations  are  given  in  Table  I. 

The  total  pressure  fluctuations  due  to  the  passage  of  a  single  vortex 
can  be  expressed  in  a  similar  manner.  Total  pressure  is  defined  as  the  sum 
of  the  static  and  dynamic  pressures.  Using  only  the  axial  component  of  velocity, 
the  total  pressure  as  measured  at  the  probe  is 

L  2  —  2  2 

PT  =  P  +  2  (U  +  u)  =  P  +  2  (U  +  2U  u  +  u  )  (5) 

1  o  o  o 

The  fluctuation  in  total  pressure,  normalized  by  velocity  head,  is  dependent  on 
the  velocity  fluctuations  and  is  given  by 


(6) 


This  total  pressure  fluctuation  is  of  interest  because  it  will  be  used  later  to 
define  the  vortex  properties  (size  and  strength)  mean  distributions. 

Statistical  Description  of  Turbulent  Flow 

Turbulent  fluid  motion  as  defined  by  Hinze,'*  is  an  irregular  condition 
of  flow  in  which  various  measurements  of  pressure,  velocity,  etc.,  show 
a  random  variation  with  time.  This  random  motion  is,  however,  amenable  to 
description  by  the  laws  of  probability  and  distinct  average  values. 

The  root  mean  square  (RMS)  level  and  power  spectral  density  (PSD)  function 
of  total  pressure  fluctuations  are  two  examples  of  average  values  that  are 
commonly  used.  Since  in  this  model  of  turbulence,  low  pressure  regions  are 
caused  by  vortices  passing  an  arbitrary  geometric  station,  it  is  necessary  to 
relate  the  RMS  level  and  PSD  function  to  the  characteristic  properties  of  the 
vortices.  These  relationships  are  developed  in  Reference  1  and  include:  (1) 
the  relationship  of  the  root  mean  square  level  to  the  mean  strength  of  the 
vortices,  and  (2)  the  relationship  of  the  power  spectral  density  function  to 
the  mean  size  of  the  vortices.  The  power  spectra  computed  from  the  analysis 
are  compared  with  that  obtained  from  test  data  to  provide  verification  of  the 
turbulent  flow  model. 
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Typical  statistical  measurements  for  a  single  probe  are  shown  in  Figure  7- 
figure  7(3)  shows  the  time  trace  of  the  total  pressure  fluctuations  which  are 
random  in  nature  -  in  the  magnitude  duration,  and  direction  of  the  pressure 
pulses.  The  time  history  can  be  expressed  as  a  probability  density  function 
(histogram)  of  the  pressure  fluctuation  as  shown  in  Figure  7(b).  Mathematically, 
the  histogram  is  the  autocorrelation  function  of  the  pressure  pulse  at  the 
probe.  The  probability  density  function  represents  both  the  range  of  pressure 
fluctuations  and  the  relative  occurrence  of  each  level.  The  time  average  is  zero, 
fiance  steady  state.  The  total  pressure  RMS  level  is  the  mean  value  of  these  pres¬ 
sure  fluctuations.  Transformation  of  this  information  from  the  time  domain  to 
the  frequency  domain  yields  the  power  spectral  density  function,  Figure  7(c) . 

These  two  functions,  Figure  7(b)  and  7(c)  define  the  vortex  properties. 

The  magnitude  of  the  pressure  fluctuation  in  Figure  7(a)  depends  directly  on  ill 
vortex  tangential  velocity.  Figure  6.  The  magnitude  of  the  fluctuations  deter':  'n 
the  histograms'  RMS  level.  Therefore,  the  RMS  level  is  a  direct  measure  ol  t  !>■ 
vortex  tangential  velocity  or  strength.  The  duration  of  a  pressure  pulse  in 
Figure  7(a)  depends  directly  on  the  size  of  the  vortex  passing  the  probe  Io.j  i 
Figure  6.  The  pressure  pulse  from  a  large  vortex  will  have  a  greater  durutioi 
than  the  pulse  from  a  small  vortex.  These  periods  of  pulse  durations  are 
expressed  as  frequency  of  occurrence  in  the  frequency  domain,  the  PSD  function. 
Figure  7(c).  The  PSD  curve  represents  the  distribution  of  vortex  sizes  tha.  p.. 
the  probe  stations.  Vortices  of  small  size  will  appear  at  the  higher  frequ.sn  : 
and  vortices  of  large  size  will  appear  at  the  lower  frequencies.  The  PSD  1  .a.  »  > 
shape  therefore,  determines  the  average  vortex  size. 

These  two  relationships  can  be  used  to  determine  the  distribution  o( 
size  and  strength  at  each  point  across  the  channel.  Unsteady  pressure  data  o 
these  forms  was  not  available  for  the  present  feasibility  study,  but  could  b 
obtained  in  a  subsequent  investigation. 
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2.2  UNSTEADY  VELOCITY  CORRELATIONS 


The  Navier-Stokes  Equations  are  considered  valid  for  turbulent  flow  if 
actual  velocities,  etc.,  are  used.  Since  the  actual  turbulent  flow  quantities 


are  functions  of  time,  the 
flow  take  the  form: 


,3U  * 

U3U 

3U 

1- —  + 

3x 

Vr—  + 

3y 

3t 

i3v  . 

3V 

Vx—  + 

3y 

5T 

Navier-Stokes  Equations 

Up-  +  u?2u) 

P  9x 

+  Up-  +  yv2v) 

P  3y 


r  two-dimensional  channel 

(7a) 

(7b) 


The  dependent  variables  can  be  replaced  by  the  mean  time  average 
fluctuation  components.  For  instance,  U  and  V  are  replaced  by  U 
The  Navier-Stokes  equations  then  take  on  the  following  form,  the 
Equations  of  Motion.  Q  _ 

a?Ux  3u2 

,  +  "  57" 

negl i gi bl e  in 
Turbulent  Flow 
\ 


3x 


+  V 


3U 


I  3P 
p  37  + 


u37  + 


3V 

'57 


]_  3P 

P  57 


=  -  —  +  v 


,3M'  3  V\ 

{Tx>'  +  3y^ 


3uv 

nr 


3v2 

w 


plus  the 
+  u  and  V 
Reynol ds 


3uv 

V* 


Reynolds  Stress  Terms 


+  v. 


whe  re  U ( t )  =  U  +  u 
V(t)  »  V  +  v 

The  upper  case  symbols  denote  the  time  average  (or  steady  state  value)  and  the 
lower  case  symbols  denote  the  fluctuation  components.  Reynolds  stress  terms  are 
dependent  on  fluctuating  components  of  the  turbulent  flow.  The  second  term  on 
the  rlghthand  side  of  the  equations  represent  the  laminar  flow  shear  stress 
terms.  Whereas,  these  terms  are  not  identically  zero,  they  are  small  compared 
to  the  turbulent  shear  stress  term  and  therefore  can  be  neglected. 


Time  averaged  unsteady  velocity  components  of  the  turbulent  flow  can  be 
established  from  the  turbulent  flow  model,  as  listed  In  Table  II.  The  axial 
and  lateral  velocity  components  are  first  described  for  a  single  vortex.  These 
are  dependent  on  time  (t)  and  on  the  vortex  properties  of  size  (a),  strength  (v  ) , 
and  spin  direction  (n) .  Unsteady  velocity  correlation  terms  are  then  found  as 
the  correlation  term  for  a  single  vortex,  multiplied  by  the  total  number  of 
vortices  per  unit  time  (N)  and  weighted  by  the  respective  probabilities  of  all 
the  random  variables.  These  are  shown  in  Table  II  in  functional  form. 


2.2.1  Solution  Procedures 

Three  methods  for  obtaining  solutions  of  the  unsteady  velocity  component 
integral  equations  in  Table  II  are  discussed.  The  first  is  the  simplest  and  is 
based  on  using  uniform  vortex  properties  across  the  channel.  The  second  method 
introduced  Is  the  most  complex.  It  was  developed  in  conjunction  with  the  use 
of  the  Turbulent  Kinetic  Energy  (TKE)  Equation.  The  final  approach  is  based  on 
using  vortex  stability  criteria  as  the  basis  for  establishing  the  vortex  flux  (N) . 
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TABLE  II 


VELOCITY  COMPONENTS  FROM  VORTEX  MODEL 


u!z'  =  uY.-t)  (tiir)  <P(Y)(P(n)(P (&*>) (P (A)  c/td* Jndv 


,n- 


(n-i) 

5  ^(fff  <?&)<?(»)  ?&*)(?(*■) 

( 1 1-2) 

nv  s  A/ jJJJ M Y) (P( *i) (Pfam ) P  (a)  dy 

( 1 1  -  3 ) 

u  »  A/  u  it)  (P  (Y )  6*  (m)  <P(fom)  (p  (&,)  dtc/O'dib/velvi  d  Y 

( 1 1  -  4 ) 

V  »  a/  v(i)  6?(y)(P{n)  (P(\Z&  *)<?(&)  dtd(Lj\/&m<inAY 

( 1 1-5) 
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2. 2. 1.1  Uniform  Properties 

The  simplest  solution  Is  obtained  by  assuming  that  the  vortex  properties 
are  constant  (uniform).  For  instance,  the  mean  velocity  profile  for  the  2-0 
channel  flow  can  be  found  from  the  integral  equation: 

*-“jG00C  u(t)P(ve)P(Y)P(a)dtdVedYda  (8) 

where,  N  Is  the  eddy  flux  (number  of  vortices  per  unit  time). 

P(  )  are  the  probability  density  functions  for  the  vortex  random  properties. 

u(t)  Is  defined  by  equation  (1-1). 

The  multiple  Integral  equation  can  be  evaluated  sequentially  by  integrating  one 
variable  at  a  time.  The  first  variable  is  the  time  variable.  The  limits  of 
integration  are  from  t  -  to  t  «  +».  That  Is,  the  Influence  of  each  eddy  Is 
considered  from  a  far  distance  upstream  to  a  far  distance  downstream  of  the 
sensing  element. 

The  second  variable  of  Integration  is  the  eddy  strength  (Vo).  This  has  been 
assumed  to  be  constant  across  the  channel.  The  minimum  strength  physically  pos¬ 
sible  Is  zero.  The  maximum  postulated  strength  In  the  eddy  model  is  the  local 
velocity  U. 

The  eddy  size  Is  also  assumed  to  be  constant.  The  minimum  size  possible  is 
zero  and  the  maximum  size  permitted  is  the  channel  height  H. 

The  eddy  location  and  spin  direction  are  treated  in  conjunction  with  one 
another.  The  eddy  location  probability  density  function  Is  assumed  to  be  uniform. 
All  the  eddies  In  each  half  of  the  channel  are  assumed  to  spin  in  the  same  direc¬ 
tion.  The  mean  velocity  at  a  probe  position  is  found  by  integrating  equation  (8) 
from  lower  wall  to  upper  wail.  Since  all  eddies  below  the  centerline  will  rotate 


clockwise,  In  this  range  the  limits  of  Integration  are  -h  to  H/2-h.  In  the  upper 
half  of  the  channel,  the  eddies  will  spin  counter-clockwise  and  the  limits  of 
integration  are  H/2-h  to  H-h.  The  resultant  integral  solution  is: 
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Equation  (9)  was  evaluated  assuming  different  values  of  the  eddy  size  a/H. 

The  velocity  profile  was  normalized  so  that  the  velocity  is  zero  at  the  wall  and 
maximum  at  centerline.  The  results  are  shown  In  Figure  8.  These  profiles  resemble 
laminar  flow  profiles  more  than  the  turbulent  profile  represented  by  the  test  data 
of  Laufer.  The  mismatch  In  profiles  is  due  to  the  uniform  distributions  assumed 
for  the  vortex  random  properties. 

2. 2. 1.2  Vortex  Started  Impulsively 

All  preceding  work  with  the  turbulent  flow  model  assumed  that  the  vortex 
properties  remained  constant  during  the  time  period  of  interest.  That  is,  the 
flow  field  about  each  vortex  was  Independent  of  time.  In  reality,  the  eddies  that 
exist  in  the  turbulent  flow  are  changing  in  size  and  strength  as  they  are  convec- 
ted  downstream.  Therefore,  the  more  basic  form  of  the  flow  field  about  a  vortex 
started  Impulsively  at  time  t=0  was  used  to  describe  the  turbulent  flow.  Each 
such  vortex  would  then  be  allowed  to  grow  In  size  and  decay  In  strength  with  time. 
However,  the  development  led  to  a  highly  intractable  integral  form. 

G.  I.  Taylor  in  1918  gave  the  solution  of  a  vortex  flow  field  which  satisfies 
the  unsteady  Navier-Stokes  Equations  and  the  proper  boundary  conditions.  This 
solution  represents  a  vortex  formed  instantaneously  In  undisturbed  flow.  The 
influence  of  this  impulse  begins  to  propogate  outward  at  time  t=0.  The  vortex 
flow  field  is  defined  by: 


r  '  X* 

ve  “  8  “2  e 

where:  Ve  =  velocity  i 

in  angular  direction 

B 

=  constant 

r 

=  radius 

t 

=  time 

V 

=  kinematic 

vi scos i ty  *  y/p 

(10) 


The  radius  (a)  at  which  the  velocity  is  maximum  can  be  determined  by  setting  the 
derivative  9V«/3r  equal  to  zero,  or: 


3v 

3r 


0 


so  that: 


r 


a 


on 


According  to  (11),  the  vortex  grows  In  size  while  according  to  (10),  it  decays 
In  strength  with  time. 
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The  vortex  tangential  velocity  can  be  expressed  as  axial  and  lateral  com¬ 
ponents  by  using  the  following  transformations: 

Vortex  Coordinates 


u  =  -v0  sin  8  =  -  ven^- 

.  x 

v  =  v#  cos  9  =  v0n— 


u 


-Bn 


2 

x  +y 


2 


**Vt 


(12a) 


v 


2  2 
x  +y 


02b) 


These  velocity  components  are  then  transformed  to  the  coordinate  system  "fixed" 
to  the  channel. 


u 


(+Xo+Ut)2+Y2 
4vt  ” 


r 


-Bn 


+Xo+Ut 
- - —  e 


(+Xo+Ut)Z+Y2 

5vt 


where:  x  ■ 
Xo  “ 


-Xo  +  Ut 

station  where  vortex  originates 


(13a) 


(13b) 


t 

t 

t 
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The  closed  solution  of  the  Reynolds  Equation  of  Motion  depends  on  definition 
of  the  cross-velocity  stress  term,  u'v‘.  Using  the  above  unsteady  velocity  com¬ 
ponent  definitions,  Equation  (13),  the  cross-velocity  correlation  term  is: 

- (Xo+Ut) 2+Y2 

u-v1  -  (Xo  +  Ut)  Y  e  2vt  (1M 

t* 

The  cross-velocity  stress  term  Is  obtained  by  integrating  with  respect  to  all  the 
vortex  properties.  The  vortex  properties  include  size,  strength  (B) ,  spin  direc¬ 
tion,  and  location  (X  and  Y) ,  and  has  the  form: 


u'v1  ■  N  fff  f  uv(t)  dtdadve  dndY 

m 


(15) 


The  x  location  Is  defined  in  terms  of  time.  Since  the  vortex  exists  only  for 
t  ^  0,  the  time  Integral  has  the  form: 


uv 


1 


Yl  ( 


Xo+Ut 


(Xo+Ut)2+Y2 

5vt 

e  dt 


(16) 


where:  T  tends  to  <*>.  This  form  Is  not  readily  integrable.  Therefore,  this 
solution  procedure  was  pursued  no  further. 


2. 2. 1.3  Constant  Properties  with  Stability  Criteria 

The  turbulent  flow  model  represents  a  totally  random  turbulent  flow  when 
the  vortex  properties  are  described  as  random  variables.  At  any  point  in  the 
flow  field  (x,  y)  the  vortices  can  have  a  range  of  size,  strength,  etc.  For 
instance,  the  vortex  size  can  vary  from  zero  to  a  channel  half-width.  The 
variation  of  such  a  random  variable  can  be  defined  by  use  of  probability  density 
function.  An  alternate  approach  for  specifying  the  vortex  properties  is  to  define 
only  the  mean  value  of  each  vortex  property  at  a  given  point  (x,  y)  in  the  flow 
field.  This  is  equivalent  to  using  only  the  mean  value  of  the  vortex  property 
rather  than  Its  probability  density  function.  It  is  necessary  to  define  only 
the  mean  value  variation  across  the  channel.  Distributions  investigated  are  pre¬ 
sented  in  Section  2.3- 


The  velocity  component  Integral  equations  in  Table  II  are  simplified  by 
using  constant  values  for  the  several  vortex  properties.  The  value  of  the  vortex 
property  is  allowed  to  vary  across  the  channel,  but  it  will  be  constant  at  any 
one  point.  That  is,  no  fluctuation  in  vortex  properties  is  allowed  at  a  given 
flow  field  point.  For  instance,  in  functional  form,  the  vortex  size  expression 
is  simplified  from  the  random  variable: 

P (a)  -  f (h/H) 


to  the  algebraic  variable: 
a  -  f (h/H) 
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The  vortex  size  therefore  can  be  defined  in  terms  of  the  distance  across  the 
channel . 

The  integral  equation  in  Table  If  begins  with  the  fluctuation  of  a  single 
vortex.  The  time  interval  represents  the  time  averaged  influence  for  a  single 
vortex.  The  vortex  spin  direction  is  chosen  as  clockwise  in  the  lower  half  part 
of  the  duct  and  counter-clockwise  in  the  upper  half,  see  Figure  6.  The  basis 
for  this  assumption  is  the  observation  that  a  velocity  decays  below  the  average 
velocity  near  the  wall,  but  Is  above  average  at  the  channel  centerline.  The 
vortex  strength  (Vem)  and  size  (a)  are  defined  as  functions  from  the  lower  wall. 

The  last  unknown  variable  is  the  vortex  flux  (N) .  The  flux  rate  can  be 
specified  by  using  a  stability  criteria,  as  shown  in  Table  lil.  The  final 
expression  (I  I  1-1)  relates  the  flux  to  the  local  velocity  (ll)  and  vortex  size 
(a)  variation  across  the  channel. 

The  development  for  the  mean  velocity  profile  is  presented  in  Table  IV, 

This  is  the  solution  of  equation  (1 1-1))  using  equation  (1-1)  for  the  perturba¬ 
tion  in  the  axial  direction.  The  flux  rate  expression  (III-!)  is  also  used. 

The  final  result: 


must  still  be  integrated  with  respect  to  distance  across  the  channel.  The  inte¬ 
gration  limits  are  shown  in  Table  I  where  Y/H  is  the  variable  of  integration. 

The  variation  of  vortex  size  (a/H)  and  strength  (Vem/Uo)  are  to  be  specified  as 
functions  of  h/H,  the  distance  from  the  lower  wall.  Typical  results  obtained 
are  discussed  in  the  next  section. 

The  integral  equations  for  three  unsteady  velocity  correlation  terms  are 
presented  in  Table  V.  These  derivations  include  the  root-mean-square  levels  of 
the  axial  (►IP2  /Uo)  and  lateral  (  ►/vT2  /Uo)  velocity  perturbations  and  the  cross 
correlation  term  u 1  v 1  /Uo2 .  The  cross  correlation  term  is  an  odd  function  of  (^i.) 
and  therefore  yields  zero.  This  assumes  that  time  averaging  the  sum  of  vortex 
influence  is  the  same  as  summing  the  vortex  influence  and  then  taking  the  time 
average.  A  finite  answer  can  be  obtained  by  considering  the  absolute  value  of 
the  function  with  time.  Typical  results  will  be  presented  in  the  Data/Analysis 
sect i on . 

Integral  equations  were  also  developed  for  the  unsteady  terms  in  the  Turbu¬ 
lent  Kinetic  Energy  equation  (Table  Vl).  This  topic  is  discussed  later  herein. 
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TABLE  III 


Number  of  Vortices  in  2-D  Duct 


Assume,  that  on  the  average,  vortices  will  align  themselves  as  shown  schemati¬ 
cally  above,  with  spacings  of  /  and  d.  Then, 

Number  of  vortices/unit  area  = 

Number  of  vortices/unit  Time  =  AN  =  *  U  Ay 


The  spacing  ratio,  l/d,  can  be  obtained  by  imposing  a  stability  criteria 
on  the  vortex  rows  as,  for  example,  in  Mi Ine-Thompson 1 s  Theoretical  Hydrodynamics, 
p.  378,  where, 


t_ 

d 


3.56 


Further  assume  that  /<*  a,  Csa.  Then  AN  =  A^-  and  in  the  limit: 

r 

...  _  3.56  U  dy 


(III-l) 


TABLE  IV 


Time  Averaqe  Velocity  Induced  By  Vortices 


Velocity  induced  by  vortex  "i"  at  probe: 

“tit)  -  *  (?l)2  -11 

The  time  averaqe  velocity  induced  will  be: 


(I-D 


VT.|  -0 


(IV-1) 


The  tine  averaged  velocity  induced  by  all  vortices  will  be  N  times  that  induced 
by  a  single  vortex,  thus: 


u  =  /  du  =  /  Uj  T  dN 


du 


3.56  u 


3.56  u  Y  pr  a  -1/3  { (-)  -H 
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Simplifying,  the  average  velocity  becomes: 
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TABLE  V 


VELOCITY  CORRELATION  TERMS 
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TABLE  V 


VELOCITY  CORRELATION  TERMS  (CONTINUED) 


USING  ACTUAL  VALUE: 
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Note:  This  assumes  that  time  averaging  the  sum  of  vortex 
influence  is  same  as  summing  the  vortex  influence  and 
taking  time  average.  This  is  not  true  and  thus  answer 
is  too  low. 
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USING  ABSOLUTE  VALUE: 
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Note:  This  assumes  influence  of  one  vortex  on  another  is 
very  (all!)  important.  Thus  answer  is  probably  too  high. 


TABLE  VI 

TURBULENT  KINETIC  ENERGY  EQUATION  TERMS 
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2.3  DISTRIBUTION  OF  VORTEX  PROPERTIES  FOR  TWO-DIMENSIONAL  FLOW 


The  turbulent  flow  node  I  Is  composed  of  a  series  of  random  vortices  con- 
vected  by  the  local  flow.  The  flow  about  each  individual  vortex  Is  essentially 
a  steady  state  flow  field  as  described  in  a  coordinate  system  fixed  to  the  vor¬ 
tex.  The  flow  field  for  each  vortex  extends  across  the  entire  channel.  The 
single  vortex  circular  motion  appears  as  a  fluctuation  or  turbulence  as  it 
passes  a  point  fixed  to  the  channel.  The  flow  at  this  observation  point  is 
influenced  by  the  flow  field  from  every  vortex  that  passes  down  the  channel. 

In  other  words,  the  turbulence  at  a  point  in  the  channel  depends  on  the  char¬ 
acteristics  of  vortices  all  across  the  channel.  Therefore,  it  is  necessary  to 
define  the  distribution  of  vortex  properties  across  the  two-dimensional  channel. 

2.J.1  Vortex  Distributions  from  Measured  Unsteady  Pressures 

As  described  in  Section  2.1,  the  Turbulent  Flow  Model  can  be  used  to  compute 
the  distribution  of  vortex  properties  Measurement  of  the  unsteady  pressure  fluc¬ 
tuations  can  yield  a  root-mean-square  (RMS)  level  and  power  spectral  density  (PSD) 
function  at  various  points  across  the  channel.  The  RMS  is  a  measure  of  the  aver¬ 
age  vortex  strength  and  the  PSD  shape  determines  the  mean  vortex  size.  This  pro¬ 
cedure  requires  measurement  of  the  turbulence  in  terms  of  the  unsteady  total  pre- 
sure  fluctuations. 

This  type  of  data  was  not  available  for  the  current  study.  Therefore,  the 
primary  focus  was  at  finding  the  typical  distributions  of  the  vortex  properties 
that  yielded  unsteady  velocity  correlation  terms  that  matched  the  test  data  of 
Laufer,  Reference  2. 

2.3.2  Distribution  of  Vortex  Properties  for  Analytical  Solution 

The  mean  velocity  profile,  Equation  (1/),  and  also  the  turbulence  terms 
u : str i but  ions  depend  on  the  distribution  of  vortex  properties  across  the  channel, 
lor  a  fully  exact  treatment,  each  vertex  property  should  be  defined  by  a  proba¬ 
bility  density  function.  This  permits  each  property  to  have  any  value  within 
:ts  acceptable  range  with  assigned  probabilities  for  each  value.  However,  for 
1  In  s  analysis,  each  vortex  property  is  assumed  to  have  only  one  value,  its  average* 
vu  I  ue ,  at  each  location  across  the  channel.  Each  vortex  parameter  is  discussed 
hr  i  e f  1  y  . 

The?  vortex  axial  location  is  assumed  to  be  uniform.  This  paramete'r  is 
1  presented  by  the  Ut/a  term  in  equation  (I-!)  and  is  included  as  the  dt  variable 
■i  integration  in  Equation  (17).  Each  eddy  is  convected  by  the  local  velocity,  so 
that  X  =  Ut . 

The  vortex  spin  direction  is  also  floated  as  being  constant  with  the  following 
i.i.nd  i  t  i  ons .  In  the  lower  half  of  the  channel  all  eddies  are  assumed  to  spin  in 
the  clockwise  direction  (n  =  -I)  and  in  the  upper  half  of  the  channel  all  vortices 
spin  counter-clockwise  (n  *»  +1);  see  Table  I. 

The  vortex  flux  rate  expression  was  derived  using  the  M i 1 ne-Thompson  sia- 
I i I i ty  criteria  which  is  based  on  rows  of  vortices  of  alternating  spin  direction. 
This  is  applicable  because  the  vortices  in  the  turbulent  flow  are  random.  The 
1  <;str iction  in  spin  direction  imposed  in  the  preceding  paragraph  is  similar  to 
using  the  average  value  of  vortex  size  and  strength. 


25 


The  vortex  maximum  tangential  velocity  (  @m/Uo)  Is  described  as  a  function 
of  location  across  the  channel  (h/H) .  Four  velocity  distributions  used  early 
In  the  study  are  shown  In  Figure  3.  The  vortex  strength  Is  expressed  as  the 
ratio  of  maximum  tangential  velocity  (V0m)  to  the  channel  maximum  velocity  (Uo) . 
Referring  to  Table  1,  the  distance  from  the  lower  wall  Is  h  and  the  channel 
width  (height)  Is  H.  The  velocity  distributions  are  shown  from  the  wall  to  the 
channel  centerline  only  and  are  assumed  to  be  symmetrical  about  the  channel 
centerline.  Distribution  VI  is  similar  to  the  classical  turbulent  boundary 
layer  shear  stress  distribution.  The  velocity  Is  zero  at  the  duct  centerline 
and  Increases  linearly  toward  the  wall,  then  decays  rapidly  to  zero  at  the  sur¬ 
face.  Distribution  V2  Is  more  or  less  arbitrary,  It  has  a  finite  value  at  the 
centerline  but  Is  zero  at  the  wall.  Distribution  V3  Is  directly  proportional 
to  the  derivative  of  the  Laufer  data  mean  velocity  profile  and  Is  not  zero  at 
the  wall,  and  \lk  Is  a  uniform  distribution.  Additional  vortex  strength  distri¬ 
butions  were  developed  during  the  study  and  will  be  presented  In  the  DATA/ANALYSIS 
COMPARISON  section  2.5. 

The  vortex  size  (a/H)  Is  also  described  as  a  function  of  location  across 
the  channel  (h/H).  Four  vortex  size  distributions  used  early  In  the  study  are 
shown  In  Figure  10.  The  vortex  size,  (If)  ,  is  presented  as  a  ratio  of  the  channel 
width  (H) .  And  the  distribution  is  specified  only  from  the  lower  wall  to  the 
channel  centerline  and  Is  also  assumed  to  be  symmetrical  about  the  centerline. 

The  A2  distribution,  with  ~S  / H  -  0.5,  Is  similar  to  that  used  previously  for  the 
closed  solution.  Distribution  A1  is  completely  arbitrary  but  reflects  the  empir¬ 
ical  notion  that  the  vortices  are  small  near  the  wall  and  Increase  in  size  toward 
the  channel  centerline.  Distribution  A3  was  derived  as  a  first  order  match  of 
the  Laufer  mean  velocity  profile  data.  Distribution  A4  is  the  average  of  A1  and 
A3.  Again,  other  vortex  size  distributions  were  developed  during  the  study  and 
discussed  in  the  DATA/ANALYS I S  COMPARISONS  section. 

The  last  variable  in  Equation  (17)  is  the  lateral  location  Y.  Note:  See 
model  In  Table  I,  that  the  variable  Y  is  measured  from  some  "sensing"  point 
located  across  the  channel.  The  range  of  values  varies  from  -h  to  H-h.  The 

parameter  h  specifies  a  point  in  the  channel  and  corresponds  to  the  location  of 

the  "sensing"  point.  The  vortex  size  (a/H)  and  strength  (V0m/U)  distributions 
are  specified  in  terms  of  h/H.  The  variable  Y  then  defines  the  contributions 
from  vortices  all  across  the  channel  at  point  h.  The  vortex  size  and  strength 

distributions  are  specified  only  for  one-half  of  the  channel  (the  other  half  is 

assumed  to  be  symmetrical).  The  lateral  position  Integration  variable  (Y)  ,  how¬ 
ever,  extends  from  wall  to  wall.  This  is  required  because  all  vortices  In  the 
channel  contribute  to  the  turbulence  at  each  point  in  the  flow  field.  This  range 
of  the  Y  variable  permits  the  vortex  spin  direction  to  change  at  the  channel 
center  1 tne. 

The  above  solution  procedure  for  the  mean  velocity  profile  is  also  applicable 
to  the  other  turbulence  terms.  Typical  results  obtained  are  shown  In  the  DATA/ 
ANALYSIS  COMPARISONS  section. 
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Figure  10.  Vortex  Size  Distributions. 
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2.1*  EFFECTS  OF  BOUNDARY  CONDITIONS 


The  objective  of  the  current  study  Is  to  obtain  a  solution  of  the  Reynolds 
Equation  of  Motion  for  turbulent  flow  using  a  statistical  vortex  turbulence 
closure.  This  link  Is  provided  by  the  Turbulent  Flow  Model,  developed  In 
Reference  1  and  applied  herein.  The  Reynolds  Equation  defines  the  steady  state 
flow  field  In  terms  of  both  steady  state  and  time  averaged  unsteady  velocity 
components.  Historically,  the  unsteady  velocity  components  have  been  related 
by  use  of  the  "mixing  lengths"  or  "eddy  viscosity".  The  Reynolds  stresses  are 
defined  In  terms  of  these  "constants".  In  the  current  study,  the  Reynolds 
stresses  are  defined  directly  In  terms  of  the  turbulent  flow  characteristics. 

The  solution  procedure  based  on  the  vortex  started  Impulsively,  Section  2.2. 
1.2,  required  an  additional  equation  in  order  to  determine  the  vortlcity  strength 
(B).  The  Turbulent  Kinetic  Energy  Equation  was  formulated  for  this  initial  use. 
Although  this  solution  procedure  was  not  completed,  solutions  for  the  TKE  terms 
were  conducted  and  are  also  discussed  in  this  section. 

2.4.1  Reynolds  Equation  of  Motion 

The  Reynolds  Equations  of  Motion  for  fully  developed  two-dimensional  flow 
are  derived  from  Section  2.2.  For  this  flow,  gradients  in  the  axial  direction  of 
the  mean  flow  and  turbulence  quantities  vanish.  Application  of  the  continuity 
equation  for  the  mean  flow  and  the  impervious  wall  boundary  condition  shows 
that  the  mean  transverse  velocity  components  are  zero  across  the  channel.  This 
removes  the  continuity  equation  from  the  set  of  governing  equations.  The 
resultant  simplified  governing  equations  may  be  written  in  dimensionless  form 
as : 


_L  ± JL L_  (iTV")  =  i  -i£ 

Re  3  y2  *  y  v  '  2  3  x 

o 


where  the  dimensionless  variables  are  defined  as  follows: 


(18) 

( I ') ) 


where,  u  and  p  are  reference  quantities  and  an  asterisk  denotes  a  dimensional 
quantity?  Equation  (19)  may  be  integrated  to  yield  the  following  algebraic 
equation  for  the  pressure  distribution  across  the  channel: 
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(21) 


•  w 


p  “  p  -2(v'2) 

W 

where  p  is  the  local  (dimensionless)  wall  static  pressure.  Differentiating 
equation  (21)  yields) 
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Substituting  this  result  in  equation  (18)  and  integrating  yields  the  following 
result  upon  application  of  the  symmetry  conditions  for  the  mean  flow  and 
fluctuating  components 
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Equat i on  (22)  may  be  integrated  (numer i ca 1 1 y )  with  the  boundary  conditions  u  = 
u 1 v‘  ■  0  at  the  wall  if  an  expression  for  tr v 1  as  f(y)  is  known.  Such  a  functional 
expression  can  be  obtained  from  the  Turbulent  Flow  Model  Solution.  Typical  results 
with  this  numerical  procedure  are  given  below. 


At  present  the  wall  pressure  gradient  must  be  known  in  order  to  integrate  the 
equation.  This  is  desirable  since  the  u ' v '  distribution  will  be  unconstrained 
during  the  turbulence  model  development.  When  a  satisfactory  model  is  obtained, 
the  pressure  gradient  term  can  be  iterated  to  meet  the  equivalent  conditions: 
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dy 
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at  t  he  center  line. 


A  computer  program  was  written  to  solve  equation  (22).  Since  the  equation 
is  simply  an  ordinary  differential  equation  for  the  mean  velocity  U,  any  suit¬ 
able  numerical  algorithm  may  be  used.  The  present  computet  program  uses  the 
multi  step  method  by  Shampine  and  Gordon,  reference  6.  The  program  was  checked 
for  a  laminar  and  turbulent  flow  case. 

The  laminar  flow  is  obtained  by  setting  u 1 v '  “  0.  An  exact  analytical 
solution  exists  for  this  flow  and  is  given  by: 
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A  comparison  of  the  analytical  and  numerical  solutions  in  Figure  11  for 
Re0  =  30800  shows  that  the  solutions  are  identical.  Although  a  laminar  flow 
probably  could  not  be  maintained  physically  at  this  high  Reynolds  number,  this 
does  not  affect  the  solution  of  the  equation. 

The  turbulent  flow  case  considered  uses  the  data  of  Laufer  for  a  Reynolds 
number  of  24600  (or  a  half  width  Reynolds  number  of  12300).  The  u 1 v 1  distri¬ 
bution  was  interpolated  by  cubic  splines  from  the  measured  data  that  was 
obtained  by  hot  wire  anemometry.  The  resultant  calculated  mean  velocity  profile 
is  compared  to  the  measured  profile  in  Figure  12  and  shows  excellent  agreement. 
Initial  attempts  to  calculate  the  mean  velocity  profile  shown  in  Figure  12 
revealed  extreme  sensitivity  to  the  u ' v 1  profile,  especially  in  the  near  wall 
region  where  gradients  are  very  large. 

The  Reynolds  stress  term  is  defined  by  equation  (l  I  - 3 ) ,  Table  II,  in  terms 
of  the  Turbulent  Flow  Model.  The  integral  solution  procedure  used  for  this  mean 
velocity  profile,  equation  (17),  is  also  applicable  to  the  Reyonlds  stress  term. 
Results  obtained  are  presented  in  Section  2.5  and  can  be  used  in  the  preceding 
technique  to  obtain  solution  of  the  Reynolds  Equations  of  Motion. 

2.4.2  Turbulent  Kinetic  Enerqy  (TKE)  Equation 

The  method  of  specifying  the  turbulent  flow  model  by  the  vortex  started 
impulsively.  Section  2. 2. 1.2,  required  additional  constraints  for  proper 
closure.  It  was  found  that  this  closure  could  be  achieved  by  using  the  turbulent 
kinetic  energy  equation.  And  although  the  solution  method  was  abandoned  due  to 
the  complexity  in  the  integration  procedure,  the  TKE  was  retained.  The  TKE  offers 
additional  means  of  verifying  the  Turbulent  Flow  Model  solution  for  turbulent  flow 
in  a  2-D  channel . 

The  TKE  has  the  following  form  for  the  2-0  channel  flow.2 
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(23) 


The  first  item  on  the  left  hand  side  is  the  turbulence  production  term. 
Turbulence  is  created  by  the  Reynolds  stresses  working  on  the  mean  velocity 
profile.  The  second  term  on  the  left  hand  side  accounts  for  diffusion  of  the 
turbulence  due  to  the  pressure  and  velocity  gradient.  These  two  parts  are 
combined  here  similarly  to  Laufers  approach.  This  term  accounts  for  the  trans¬ 
fer  of  turbulence  into  (and/or  from)  a  region  in  the  flow.  The  right  hand 
term  accounts  for  the  dissipation  of  turbulence.  The  vortices  that  constitute 
the  turbulence  yield  their  energy  as  they  disappear. 


The  three  terms  in  the  TKE,  equation  (23),  represent  energy  levels.  Laufer 
presents  these  terms  as  the  change  in  energy  per  unit  time  transferring  into  and 
out  from  a  unit  volume  at  each  point  across  the  channel.  A  similar  approach 
has  been  followed  herein.  The  Turbulent  Flow  Model  was  used  to  generate  the 
correlations  in  the  diffusive  and  dissipation  terms.  Results  obtained  are 
presented  in  the  OATA/ANALYS I S  COMPARISONS  section.  Since  this  energy  relation 
must  also  be  satisfied,  it  can  also  be  used  as  a  further  check  on  the  validity 
of  the  Turbulent  Flow  Model. 
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2.5  DATA/ANALYS I S  COMPARISONS 


Results  obtained  with  the  Turbulent  Flow  Model  are  compared  with  the 
measured  results  published  by  Laufer.2  First  a  brief  review  of  the  test 
and  the  type  of  measurements  made  is  presented.  Then  comparisons  are  made 
for  mean  velocity  profiles,  Reynolds  stresses,  and  TKE  terms. 

2.5.1  Laufer' s  Measurements 
2 

Laufer  used  both  a  1-Inch  and  5-inch  channel  for  measuring  the 
unsteady  velocity  fluctuations.  For  the  5-inch  tunnel  configuration  the 
recorded  measurements  at  three  Reynolds  numbers:  12,300,  30,800,  and  61,600. 

These  Reynolds  numbers  are  based  on  the  maximum  duct  velocity  (Uo)  and  the  duct 
half-width  (2.5  Inch).  A  summary  of  the  Laufer  data  Is  presented  in  Figure  13- 
In  these  figures  the  velocities  are  normalized  by  the  maximum  duct  velocity  (Uo) 
and  are  plotted  versus  distance  from  the  lower  wall  (h/H) .  The  upper  left  chart 
shows  the  measured  (symbols)  mean  velocity  profiles.  The  exact  laminar  pro  Me 
is  included  for  reference.  The  upper  right  hand  chart  shows  the  cross-correlation 
term  -  both  the  measured  values  (symbols)  and  the  values  computed  (lines)  to 
satisfy  the  mean  profiles.  The  lower  charts  show  the  unsteady  velocity  correla¬ 
tion  measured  root-mean-square  values  for  the  axial  and  lateral  directions. 

Laufer  also  computed  t he  energy  levels  for  the  production,  diffusive,  and 
dissipation  terms  of  the  Turbulent  Kinetic  Energy  (TKE)  Equation.  These  energy 
levels  are  presented  in  Figure  1 T*  for  the  Reynolds  number  =  30,800  case.  The 
production  term  Is  zero  at  the  channel  centerline  and  increases  as  the  wall  is 
approached.  The  diffusive  term  is  negative  at  the  centerline  and  also  increases 
as  the  wall  is  approached.  The  negative  sign  at  the  centerline  Indicates  a  gain 
in  energy  in  this  region.  The  dissipation  term  remains  positive  throughout  and 
decreases  toward  the  centerline.  At  any  point  in  the  channel,  the  total  energy 
must  be  conserved.  Therefore,  the  sum  of  production  less  diffusive  less  dissi¬ 
pation  is  zero,  or: 

Production  -  Diffusive  -Dissipation  =  0 

2.5- 2  Analytical  Predictions 

Some  typical  results  obtained  in  the  study  are  presented  next.  These  com¬ 
parisons  include  the  various  parameters  for  which  test  data  were  contained  in 
Laufer's  report.  Among  these  parameters  are:  the  mean  velocity  profile,  the 
Reynolds  stress  term  (u 1 v 1 )  ,  the  unsteady  velocity  correlation  terms  (  JJ'2  and 
V'2)  and  the  Turbulent  Kinetic  Energy  equation  energy  balance  terms. 

2.5- 2. 1  Mean  Velocity  Profiles 

The  velocity  perturbation  due  to  a  single  vortex  Is  given  by  equation  (l-l). 
The  turbulence  properties  depend  directly  on  the  number  of  vortices  influencing 
the  sensing  probe.  The  vortex  flux  in  a  2-D  channel  is  derived  In  Table  ill,  and 
the  time  average  velocity  induced  by  all  the  vortices  Is  derived  in  Table  IV. 

This  approach  Is  based  on  finding  the  velocity  induced  by  a  single  vortex  at  a 
point  (probe)  and  then  finding  the  average  velocity  Induced  due  to  all  the  vor¬ 
tices  contained  in  the  flux  developed  in  Table  III.  The  resultant  average  velo¬ 
city  Is  then  given  by : 
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Mean  Velocity  Profile 


Reynolds  Stress  Distribution 
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Figure  13.  Summary  of  Turbulent  Flow  Data  in  2-D  Channel,  Reference  (2) 
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Figure  14.  Turbulent  Energy  Balance  Data.  Reference  (2). 
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Equation  (17)  was  solved  by  a  numerical  Integration  procedure.  The  variations 
of  vortex  strength  (V0m/Uo)  and  size  (a/H)  were  specified  as  per  Figures  9  and 
10.  Typical  results  are  shown  In  Figure  15.  Most  of  these  resultant  mean 
velocity  profiles  resemble  the  laminar  profile.  The  A3 V I  solution  (l.e.,  by 
using  vortex  size  distribution  A3,  Figure  10,  and  vortex  strength  distribution 
VI,  Figure  9,  Indicates  that  there  are  distributions  of  vortex  properties  which 
will  yield  correct  results. 

2. 5.2. 2  Unsteady  Velocity  Correlation  Terms 

The  Reynolds  stress  term  and  the  unsteady  velocity  correlation  terms  are 
computed  from  expressions  similar  to  Equation  (17),  as  presented  in  Table  V. 

2. 5. 2. 3  Turbulent  Kinetic  Energy  Terms 

The  turbulent  kinetic  energy  equation  for  fully  developed  two-dimensional 
channel  flow  Is  given  as  Equation  (23).  The  first  term,  the  production  term, 
is  based  on  the  Reynolds  stresses.  The  second  term  is  the  derivative  of  the 
pressure  and  velocity  transport.  The  integral  expressions  for  these  terms  are 
listed  in  Table  VI.  The  right  hand  term  of  Equation  (23)  Is  the  dissipation  term 
and  is  also  specified  in  Table  VI. 

2.5  3  Comparisons 

The  combination  A2V*»  represent  the  solution  for  constant  vortex  size  and 
constant  distributions.  Results  obtained  with  these  distributions  of  vortex 
properties  are  shown  In  Figure  16.  The  mean  velocity  profile  predicted  Is 
similar  to  a  laminar  profile.  Note  that  although  the  Reynolds  stress  term  dis¬ 
tribution  is  similar  In  shape  to  the  test  data,  the  values  are  more  than  an 
order  of  magnitude  high.  The  unsteady  velocity  correlation  terms  are  also  an  order 
of  magnitude  high.  The  energy  production  and  diffusive  terms  are  three  orders 
of  magnitude  too  large.  This  Is  due  to  the  laminar  type  profile. 

The  combination  A3V1  solution  yields  mean  velocity  profile  which  Is  closer 
to  the  turbulent  profile.  The  Reynolds  stress  and  unsteady  velocity  correlation 
terms  are  reduced  to  less  than  a  factor  of  ten  high.  The  TKE  production  and 
diffusive  terms  are  reduced  to  about  an  order  of  magnitude  difference,  too. 

These  comparisons  are  shown  in  Figure  17. 

A  large  number  of  vortex  strength  and  size  distributions  were  Investigated 
during  the  study.  Six  strength  and  three  size  distributions  used  are  shown  in 
Figure  18.  The  velocity  distributions  shown  all  have  the  general  properties  of 
distribution  VI,  Figure  9.  The  effects  of  each  distribution  are  discussed  below. 
The  three  vortex  size  distributions  shown  in  Figure  18(b)  have  the  general  char¬ 
acteristics  of  distribution  A3,  Figure  10.  The  initial  slope  (dS/dh)  is  the 
principal  variable.  The  net  effect  is  to  make  A 1 50  distribution  almost  linear 
and  the  A300  distribution  full  and  flat.  The  following  trends  were  observed. 

The  magnitude  of  the  Reynolds  stress  and  the  unsteady  velocity  terms  depend 
directly  on  the  vortex  strength.  The  level  Is  more  significant  than  the  dis¬ 
tribution;  the  distribution  is  most  important  near  the  wall.  The  vortex  strength 
distribution  seems  to  have  minor  effect  on  the  mean  velocity  profile. 
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Figure  16.  Comparisons  of  Data  and  Predictions  with  Distributions  A2V4. 
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The  vortex  size  distribution  affects  the  mean  velocity  profile  and  the 
Reynolds  stress  term.  The  mean  velocity  profile  Is  sensitive  to  the  Initial 
slope  (at  the  wall)  of  the  vortex  size  distribution.  This  Initial  slope  also 
affects  the  location  of  the  peak  value  of  Reynolds  stress  term.  The  level  of 
the  Reynolds  stress  term  depends  on  the  vortex  size  Initial  value  (at  the  wall). 

These  observations  are  generally  Illustrated  by  the  following  comparisons. 
Figures  19  through  2k.  The  predictions  are  compared  to  the  Re  ■  30,800  flow 
case  test  data.  Each  figure  corresponds  to  one  of  the  vortex  strength  distri¬ 
butions  of  Figure  18(a)  and  Includes  the  three  vortex  size  distributions  of 
Figure  18(b) . 

Mean  Velocity  Profiles 

Note  that  the  Initial  slope  of  the  size  distribution  has  a  definite  effect 
on  the  profile  shape.  The  lower  the  Initial  slope,  the  fuller  the  profile. 

Vortex  strength  distribution  V3^,  V36  and  \lk 2  yield  the  best  match  In  terms  of 
the  mean  velocity  profile.  The  other  vortex  strength  distributions  yield  pro¬ 
files  that  are  too  full.  This  Is  characteristic  of  a  higher  Reynolds  number 
flow;  see  Figure  13. 

Reynolds  Stress  Term 

The  variation  of  Reynolds  stresses  computed  using  vortex  strength  distri¬ 
butions  V3i»,  V36  and  \Jk2  agree  well  with  the  data  near  the  channel  centerline 
but  deviate  the  most  towards  the  wall.  Note  that  the  peak  values  occur  further 
from  the  wall.  The  initial  slope  of  the  vortex  size  distribution  has  a  secon¬ 
dary  but  definite  effect  on  the  magnitude  of  the  Reynolds  stress. 

Unsteady  Velocity  Correlations 

The  model  now  predicts  the  unsteady  velocity  correlation  terms  (u ' 2  and 
v'*)  that  are  approximately  equal  to  the  test  data.  For  the  axial  component 
(u'*)  ,  the  predicted  slope  Is  generally  lower  than  the  test  data.  The  vortex 
size  variation  with  the  lowest  initial  slope  (A150)  generally  yields  the  best 
agreement.  For  the  lateral  component  (v 1 *) ,  the  value  near  the  centerline  is 
well  predicted,  but  tends  to  be  high  near  the  wall. 

Turbulent  Kinetic  Energy  (TKE)  Terms 

Vortex  distributions  V34,  V36,  and  \lk 2  yield  the  best  comparisons  for  the 
production  and  dissipation  terms.  The  diffusive  term  tends  to  be  a  little  high. 
For  the  other  vortex  strength  distributions,  the  production  term  tends  to  be 
slightly  lower.  In  all  these  comparisons,  the  TKE  terms  Increase  rapidly  as 
the  wall  is  approached.  This  is  due  to  the  high  turbulence  levels  near  the  wall. 
Laufer  presented  TKE  levels  only  in  the  region  away  from  the  wall. 

Summary 

Vortex  velocity  distributions  V36  and  Vk2  give  the  best  correlations  with 
the  test  data.  The  mean  velocity  profile  is  well  represented.  The  characteristic 
shape  of  the  unsteady  velocity  fluctuations  are  also  well  characterized,  as  are 
the  turbulent  kinetic  energy  equation  terms.  The  Reynolds  stress  distribution 
also  correlates  well  near  the  duct  centerline,  but  tends  to  deviate  near  the  wall. 
Overall,  the  V29B/A300  distribution  best  matches  the  total  data  set. 

The  statistical  vortex  model  of  turbulence  Is  based  on  the  hypothesis  that 
turbulence  is  composed  of  a  series  of  vortices  (eddies)  with  random  properties 
which  are  convected  downstream  by  the  mean  flow.  The  unsteady  velocity  cor¬ 
relations  (and  Reynolds  stress  term)  are  expressed  In  terms  of  the  turbulence 
model  vortex  properties.  These  properties  can  be  determined  from  measurements 


Mean  Velocity  Profile 


f/S*  *150  _ 
TT  A225 
F  A300  I 


Reynolds  Stress 

V'AflO 
A225 
A3  00 


■■■ 


—  Cow 

by 


Computed 
by  Laufer 


Velocity  Correlations 


T| 

_ a 

10! 

■ 

■ 

■ 

■1 

B 

B 

■ 

_ A30 

IHHNS1 

j 

*  V^-^JA225 
4 - ? — — - T  A1 50 


Oeta  froa  Laufer,  Reference  (2) 
Re  •  30,600 


Test  Data 
Prediction 


Turbulent  Kinetic  Energy 
Equation  Terms 

PRODUCTION 


DISSIPATION 


DIFFUSIVE 


0  0.1  0.2  0.3  0.4  0.5 

h/H 


3  0.1  0.2  0.3  0.4  0.5 

h/H 


Figure  19.  Comparisons  of  Oata  and  Predictions  with  Distribution  V298. 
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of  the  flow  field  unsteady  pressure  fluctuations.  Since  adequate  measurements 
were  unavailable,  the  vortex  properties  distributions  were  assumed  and  refined 
during  the  data/analys Is  comparisons  to  achieve  agreement  with  the  measure 
turbulence  parameters. 

Analytically,  the  statistical  model  of  turbulence  forms  the  closure  required 
for  the  governing  Reynolds  and  turbulent  kinetic  energy  flow  equations.  The 
solutions  obtained  completely  characterize  the  flow  with  a  single  distributive 
set  of  vortex  (eddy)  properties.  Comparisons  with  a  2-D  channel  data  base 
verified  the  character Izat Ions. 

2.6  ANALYSIS  FOR  AX  I  SYMMETRIC  FLOW  AND  JET  MIXING 

An  outline  Is  presented  herein  for  extending  the  Turbulent  Flow  Model  to 
the  case  of  Axisymmetric  Flow  and  Jet  Mixing 

2.6.1  Axisymmetric  Flow 

The  flow  field  for  fully  developed  turbulent  flow  in  a  duct  of  circular 
cross  section  Is  customarily  described  in  cylindrical  coordinates.  The  Reynolds 
Equation  of  Motion  for  this  case  is  well  established.  The  Turbulent  Flow  Model 
relations  typified  for  the  two-dimensional  case  are  equally- appl 1  cable  to  the 
axisymmetric  flow  case.  The  distributions  of  the  vortex  properties  would  be 
different,  but  these  distributions  can  be  established  from  measurements  and/or 
comparisons  with  available  data. 

2.6.2  Jet  Mixing 

The  j'et  mixing  case  is  fundamentally  a  free  shear  flow  field.  This  type  of 
flow  can  be  two-dimensional,  axisymmetric,  or  of  arbitrary  shape.  The  Reynolds 
Equations  of  Motion  for  the  2-D  and  axfsymmetrfc  free  shear  flow  are  also  well 
established.  The  Turbulent  Flow  Model  must  be  adapted  to  this  type  of  flow. 

2.6.3  Outline  of  Procedure 

The  following  steps  are  necessary  to  extend  the  procedure  to  axisymmetric 
flow  and  Jet  mixing: 

For  Axisymmetric  Flow  Case: 

o  Assume  axisymmetric  flow 
o  Use  cylindrical  coordinate  system 
o  Simplify  Reynolds  Equations  of  Motion 
o  Derive  Turbulent  Flow  Model  terms  with 

-  origin  at  centerline 

-  zero  velocity  at  wall/maximum  at  centerline 

-  establish  vortex  property  distributions 

-  compare  results  to  test  data 

For  Jet  Mixing: 

o  Assume  axisymmetric  flow 
o  Use  Reynolds  Equations  derived  above 
o  Adapt  Turbulent  Flow  Model  terms  for 

-  maximum  velocity  core  flow  region 

-  changing  (expanding)  boundary 

-  freestream  velocity  beyond  boundary 

-  verify  with  test  data 


3.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  feasibility  study  to  solve  the  Reynolds  Equations  of  Motion  for  fully 
developed  turbulent  flow  In  a  two-dimensional  channel  has  led  to  the  following 
conclusions: 

-  The  turbulent  flow  model  provides  a  means  for  obtaining  closure  for  the 
governing  Reynolds  Equations  of  Motion  and  Turbulent  Kinetic  Energy  Equation. 

-  The  turbulent  flow  model  yields  a  direct  solution  of  the  mean  velocity 
prof I le. 

-  The  turbulent  flow  model  yields  the  distribution  of  the  Reynolds 
stress  across  the  channel  to  solve  the  Reynolds  Equations. 

-  The  turbulent  flow  model  also  yields  good  agreement  In  the  axial  and 
lateral  unsteady  velocity  correlation  terms. 

-  The  turbulent  flow  model  yields  good  agreement  for  the  production, 
dissipation,  and  diffusive  terms  of  the  turbulent  kinetic  energy 
equation. 

The  agreement  in  all  these  aspects  of  the  turbulent  flow  using  a  single 
distribution  set  of  vortex  or  eddy  functions  confirm  the  applicability  of 
the  Turbulent  Flow  Model  to  the  solution  of  the  general  turbulent  shear  flow 
problem.  Therefore,  it  Is  highly  recommended  that  the  overall  program  of  Jet 
Scaling  be  continued.  The  next  phase  is  to  extend  the  analytical  techniques 
to  the  case  of  axisymmetric  flow  and  jet  flow. 
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NOMENCLATURE 


Al ,  A2,  etc. 

a 

B 

cs 

d 


-  vortex  size  distribution  labels 

-  vortex  core  size,  see  Figure  6 

-  vortex  model  vorticity  constant,  Equation  (14) 

-  vortex  flux  spacing  constant,  see  Table  ill,  c$ 

-  vortex  lateral  spacing,  see  Table  III 


-  2.71828 


1.845 


f<  ) 

H 

Hz 

l 

N 


P,P 
PSD 
P(  ) 
q 
q 

r 

Re ,  Reo 
RMS 
SF 
t 

TKE 
U,  V 
u,  v 


-  function  of 

-  two-dimensional  channel  height 

-  Hertz 

-  vortex  axial  spacing,  see  Table  III 

-  vortex  flux  rate,  see  Table  III 

-  vortex  spin  direction,  see  Table  I 

-  pressure 

-  power  spectral  density  function 

-  probabi 1 1 ty  of  (  ) 

-  dynamic  pressure 

2  2  2 

-  resultant  velocity,  q  =  u  +  v 

-  vortex  radius,  r  ■  a  at  maximum  velocity 

-  Reynolds  number 

-  root-mean-square  value 

-  scale  factor,  see  Figures  16  and  17 

-  time 

-  turbulent  kinetic  energy 

-  absolute  velocity  In  axial  (X)  and  lateral  (Y)  direction 

-  unsteady  velocity  components  In  axial  (X)  and  lateral  (Y)  direction 
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NOMENCLATURE  (Continued) 


u'2,  v'2 

-  unsteady  velocity  correlation  terms 

u*v' 

-  Reynolds  stress  term 

VI,  V2 ,  etc. 

-  vortex  strength  distribution  labels 

Ve 

-  vortex  tangential  velocity 

X,  Y 

-  axial  and  lateral  direction  coordinates 

Greek  Symbols 

A 

-  Increment  or  change 

e 

-  vortex  tangential  direction 

u 

-  absolute  viscosity 

V 

-  kinematic  viscosity 

ir 

-  3.14159 

P 

-  density 

Subscripts 

i 

-  individual  vortex 

m 

-  maximum  value 

o 

-  based  on  channel  maximum  velocity 

w 

-  at  the  wal 1 

e 

-  In  vortex  tangential  direction 
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